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"p^ , We present an explicit and simple form of the renormalization group equation 

Qh| which governs the quantum evolution of the effective theory for the Color Glass 

Condensate (CGC). This is a functional Fokker-Planck equation for the proba- 
bility density of the color field which describes the CGC in the covariant gauge. 
It is equivalent to the Euclidean time evolution equation for a second quantized 
current-current Hamiltonian in two spatial dimensions. The quantum corrections 
^ . are included in the leading log approximation, but the equation is fully non-linear 

with respect to the generally strong background field. In the weak field limit, it 
reduces to the BFKL equation, while in the general non-linear case it generates the 
evolution equations for Wilson-line operators previously derived by Balitsky and 
Kovchegov within perturbative QCD. 
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The problem of the gross properties of hadron interactions in the high energy limit has 
been a central problem of particle and nuclear physics for over fifty years. Remarkable 
progress in recent years has resulted from the observation that the gluon density of the 
small x part of the hadron wavefunction, the relevant piece for high energies, grows rapidly 
as x decreases |lj, |2| . When the gluon density increases, important non-linear phenomena 
are expected, which should eventually lead to saturation 0, f|, |5], 0, that is, to a limitation 
of the maximum gluon phase-space density. At saturation, this density is expected to be 
of order l/a s || 0, §, since interactions among the gluons will cutoff further growth 
once the interaction energy is comparable to the kinetic energy. For such a large density, 
perturbation theory breaks down even if the coupling constant is small, because of strong 
non-linear effects. 

On dimensional grounds, the square of the typical momentum scale associated with 
saturation is expected to be of the order of the density of gluons per unit area |J. If 
this saturation momentum Q s is large enough (Q s 3> Aq C d), then the saturation regime 
is characterized by weak coupling a s {Q 2 s ) <C 1 and large occupation numbers ~ l/a s 
which, by the correspondence principle, is a classical regime. This motivated McLerran 
and Venugopalan to describe the gluon component of the hadronic wavefunction at small 
x as a stochastic classical Yang-Mills field ||, subsequently interpreted as a "Color Glass 
Condensate" (CGC) ||. This has the advantage that the non-linear physics at saturation 
can be studied in the simpler setting of a classical field theory, via a combination of 
analytic ]5], [10], [IT], ^ || and numerical methods [0. 



In the McLerran- Venugopalan (MV) model, all the dynamical information about the 
"fast" partons (i.e., the partons with longitudinal momenta p + larger then the scale A + = 
xP + of interest^) is encoded in a probability distribution for the classical color field 
describing the "soft" gluons (i.e., the gluons with momenta k + < A + ). In principle, 
this probability distribution can be computed in perturbation theory, by integrating out 
the fast quantum modes in layers of p + , to leading order in a s ln(l/x), but to all orders 
in the strong classical background field [0. This results in a functional renormalization 
group equation (RGE) which describes the evolution of the probability density with the 
separation scale A + . The formal structure of this equation has been first presented in 
Refs. and further discussed in Refs. |14|, but previous attempts [14]] to compute 



the coefficients in this equation have suffered from technical complications and, moreover, 
appear to be inconsistent [plj with the results obtained from perturbative QCD |16|, [17]. 



In a parallel development due to Balitsky JT6[ and Kovchegov [|TjJ, the evolution of 



Wilson line operators has been studied within the operator product expansion and the 



dipole model of Mueller [|II|] . This has led to a set of coupled evolution equations for the 
expectation values of these operators |jIE| , which decouple only in the large A" limit [|Tj]] 
with A" the number of colors. (See also Refs. |19[ for some recent attempts towards solving 



5 Throughout, we consider the hadron in its infinite momentum frame, and use light-cone vector 
notations, = («+, v~, v_l), with v+ = {1/V2)(v° + v 3 ), v~ = {1/V2){v - v 3 ), and v_l = (v\v 2 ). The 
hadron four-momentum reads = (P + , 0, Oj_), with P + — > oo, while soft gluons have k + = xP + , with 
x < 1. 
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these equations.) In a recent paper P0 |, Weigert has shown that Balitsky's equations can 
be summarized as a functional evolution equation for the generating functional of the 
Wilson line operators. As we shall show in this paper, this is precisely the equation which 
governs the evolution of the probability distribution for the CGC. 

Specifically, our main purpose in this Letter is to present explicit expressions for the 



coefficients in the RGE describing the quantum evolution of the CGC [13|. The details 



of our calculation, which relies heavily on the formalism developed in Ref. ||, will be 



presented somewhere else [21 



We start by briefly reviewing the MV model and its quantum evolution (see [Q for a 
recent comprehensive discussion). In the infinite- momentum frame, and to the order of 
interest, the fast (p + > A + ) partons can be replaced by a classical color source, with color 
charge density p a {x), moving along the light-cone (z = t or x~ = 0) in the positive x + 
direction. Via the Yang-Mills equations: 

(D v F»%(x) = 6» + Pa (x) , (1) 

this source determines the classical color field A^ which describes the soft (k + ^ A + ) 
gluons. In the saturation regime, the source and the fields are strong, A ~ p ~ 1/g, and 
the classical problem is fully non-linear. 

The structure of p a (x) is largerly determined by kinematics: Since generated by shorth- 
wavelength modes, with momenta p + > A + , p is localized near x~ = 0, within a distance 
Ax~ ^ 1/A + ; it is static (i.e., independent of the LC time x + ) since the soft gluons have 
a comparatively short lifetime Ax + ~ l/k~ oc k + oc x, during which the dynamics of the 
"fast" partons can be ignored; it is a random quantity, since this is the instantaneous color 
charge of the fast partons "seen" by the shortlived soft gluons at some arbitrary time. 
The spatial correlators of the classical source p a (x) (with x = (a;~,xj_)) are inherited 
from the (generally time-dependent) quantum correlations of the fast gluons, via the 
renormalization group equation to be presented below. They are encoded into a weight 
function for p a (x), denoted as Wa[p], which is assumed to be gauge- invariant, and which 
depends upon the separation scale A + (see below). 

Physical quantities like the gluon distribution function are most directly related to 
correlations functions of the transverse fields A % in the light-cone (LC) gauge A+ = 0. In 
the MV model, the correlations functions at the scale A + are obtained as: 

(A i a (x+, x)Al(x\y) • • -) A = fvp W A [p] A'MMy) (2) 



where A\ = A\[p] is the solution to eq. ([D in the LC gauge, and is implicitely given by 
the following equations (the other components of the solution vanish: A + = A~ = 0): 

9 

U\x~,x ± ) = Pexp^igJ X dz' a a {z~, x ± )T a | , 
-Vla(x) = U\x)p{x)U{x) = p{x). (3) 
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There few important remarks about this classical field: 

a) The equations (0])-(@) above are those for a glass (here, a color glass): There is a 
source, and the source is averaged over. This is entirely analogous to what is done for 



spin glasses when one averages over background magnetic fields |22 . 

b) The solution A (x) is a two-dimensional pure-gauge (in the sense that JP- 7 = 0), but not 
also a four-dimensional pure-gauge, since the associated electric field = d + A l is non- 
vanishing. In fact, A{x) is the gauge transform with gauge function U(x) of the solution 
A 1 in the covariant gauge d^A^ = (COV-gauge), which reads simply A 1 = 5^ + a with 
a(x) satisfying: 

- Vi«(f) = p(x) , (4) 
where p(x) is the color source in the COV-gauge. 

c) Eq. (|3|) expresses the LC-gauge solution A as an explicit functional of the COV-gauge 
source p: A a = -44[/5]. On the other hand, its dependence upon the LC-gauge source p 
is not known explicitly. This, together with the gauge-invariance of the weight function 
Wa[p], suggests that it is more convenient to perform the functional averaging in eq. @ 
by integrating over p (rather than p). 

d) In writing eq. (|3]), we have chosen "retarded" boundary conditions in x~ : A l (x) — > 
for x~ — > — oo . 

e) We shall see later that the source p (and therefore also the "Coulomb" field a) has 
support at ^ x~ ^ 1/A + . Thus, when probed by fields with a poor longitudinal 
resolution (i.e., which carry momenta k + <C A + ), the configuration (0) appears effectively 
as a ^-function in x~ : 

A(x-,x ± ) w 6(x-)-V(&Vi) = 6(x-)J\L(x ± ), (5) 
9 

with A QO (x±) the asymptotic field, and 

V^(x±) = Pexp jig J dx~ a(x~, , (6) 

where the integral over x~ runs effectively up to x~ ax ~ 1/A + (since a vanishes for larger 
values of x~). On the same resolution scale, the electric field appears as a 5-function: 

F +l (x) = 8 + A w 5(a;-)^(a; ± ). (7) 



In the original MV model, eqs. (0)-(@) are simply postulated, and the weight function 
Wa[p] is not computed from first principles. Indeed, the physical origin of the large gluon 
densities, namely the quantum evolution towards small x, is absent from the model. But in 
Refs. [0, [14], one has shown that the MV model is indeed consistent with the quantum 
evolution, which entails a renormalization group equation (RGE) for the running of WaIp] 
with A. This is the same quantum evolution which in the weak field regime leads to the 
BFKL equation. But for small enough x, when one is close to saturation, this evolution 
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takes place in the strong background field A % generated at the previous steps. When 
A 1 ~ 1/g, the mean field effects are non-perturbative, and must be included exactly. 

To describe the quantum evolution, consider a sequence of two classical effective the- 
ories ("Theory I" and "Theory II") valid at the scales A + and bA + , respectively, with 
b < 1 and such as a s ln(l/&) < 1. To the order of interest, the gluon correlations at the 
lower scale bA + can be computed in two ways: As classical correlations within Theory 
II (which is valid at this scale), or by allowing for quantum fluctuations in Theory I and 
integrating out the "semi-fast" quantum gluons with longitudinal momenta in the strip 

6A+ < \p + \ < A+, (8) 

to leading order in a s ln(l/6), but to all orders in the background fields and sources. 



(See Refs. |9| for the quantum generalization of the MV model.) The difference 
AW = Wb\ — W\, and therefore the evolution equation for W\, can be then obtained 
by matching these two calculations ||. The explicit matching performed in Ref. || 
shows that it is indeed possible to absorb the effects of the semi-fast gluons on the soft 
correlation functions into a renormalization of the weight function for the effective theory. 
This is in agreement with the original conjecture in Refs. [7| and justifies the use of 
the McLerran-Venugopalan model as an effective theory at small x. 

The RGE is most naturally written in the LC-gauge, where the quantum effective 
theory is defined. It readsQ (with r = ln(P+/ A+ ) and dr = ln(l/&)) [0, [J : 

dr \2 5p*(x ± )5p b T (y ± ) 1 AxyS 5p a T (x ± ) [ xS J 

with the coefficients cr" = <7 a (x±) and Xxy = Xab(x±, y±) related to the 1-point and 2-point 
functions of the color charge Sp a (x) of the semi-fast gluons via the following relations: 

a s \n^a a (x 1 _) = J dx~ (6p a (x)) , 

a s \n^-Xab(x±,y±) = J dx~ J dy~ (5 p a (x + , x) 5 p b (x + , y)) , (10) 

where the brackets denote the average over quantum fluctuations in the background of 
the tree-level color fields A 1 . 

Note that a and x are two-dimensional densities in the transverse plane; the longitudi- 
nal structure of the relevant charge correlators has been integrated over. This is consistent 
with our picture of the quantum evolution where the semi-fast gluons are unable to dis- 
criminate (because of their low p + momenta) the longitudinal structure of the source. For 
this picture to be validated, however, one still has to check that a and x are n °t sensitive 
to the internal structure of the source (see below). 

6 We use r to denote the rapidity variable, r = ln(l/a;), rather than the more common notation y to 
avoid confusion with the space-time coordinate y. Also, we use compact notations, where repeated color 
indices (coordinate variables) are understood to be summed (integrated) over. 
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There is nevertheless a trace of this longitudinal structure in eq. (|9]): the functional 
derivatives there are to be taken with respect to the color source Pr(x±) = p a {x~ = 
1/A + , x±) in the highest bin of rapidity (r, r+dr), that is, at x~ in the range (1/A + , 1/6A + ). 
This is so since, as we shall shortly discover, the quantum corrections have support in 
that range. This longitudinal structure is important to make the functional derivatives in 
eq. (H) well defined (see also eq. ( p5|) below). 

Eq. @ is a functional Fokker-Planck equation describing diffusion in the functional 
space spanned by p with "drift velocity" a s a and "diffusion constant" a s \- Since a and x 
are generally non-linear functionals of p (see below), eq. @ is equivalent with an infinite 
set of coupled evolution equations for the correlators of p. For these equations to be 
useful, the functionals a[p] and x[p] have to be known explicitly. To this aim, it is more 
convenient to use the COV-gauge source p (rather than the LC-gauge source p) as the 
independent variable. Indeed, a and x depend upon the color source via the classical field 
A 1 which is known explicitly only in terms of p (cf. eq. (|3|)). 

As explained in Ref. 0, when going from the LC-gauge to the COV-gauge, the 
induced color charge a acquires a new contribution, in addition to the color rotation with 
matrix W . This is so because the transformation between the two gauges depends upon 
the charge content in the problem: this was p at the initial scale A + , thus giving a gauge 
transformation t/^[p], but it becomes p + 5p, with fluctuating 5p, at the new scale bA + , 
thus inducing a fluctuating component in the corresponding gauge function W [p + 5p] . 
After averaging out the quantum fluctuations, one is left with a RGE for W r [p] which is 
formally similar to eq. (||), but with p — > p, x — ► X, and a — > b, where: 

Xab(x±,y±) = Vj c (x±)xcd(x ± ,y 1 _)V db (y 1 _), b a (x±) = V^(x x ) a b (x±) - 5a a (x±), 
Sa a (x ± ) = 9 -f abc J 'd 2 y ± Xcb(x±,y±)(y±\^\x 1 _) . (11) 

The correction —5a, to be referred to as the "classical polarization" , is the result of the 
quantum evolution of the gauge transformation to the COV-gauge. 

It is our purpose in this Letter to present explicit expressions for the quantities x an d 
b as functionals of p, and thus completely specify the RGE for W T [p]. The Feynman rules 
for the calculation of these quantities have been already presented in Refs. [13, They 
involve the propagator 

iG%{x,y)[A] = (T<(zK(y)) (12) 
of the semi-fast gluons in the presence of the background fields A 1 and in the LC- 



gauge a + = (where the quantum theory is a priori formulated |T3|, || ) . In Ref. || , this 
propagator has been constructed to all orders in the background field, by exploiting the 
special geometry of the latter (see also Refs. [f23], [24]]). One of the subtle points in this 



construction has been the choice of a gauge condition (i.e., of an ie prescription for the 
pole in l/p + in the gluon propagator) to fix the residual gauge freedom in the LC-gauge. 
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This is important since the general results for x an d cr turn out to be dependent upon 
this prescript iorf]. 

For consistency with the retarded boundary conditions imposed on the classical solu- 
tion (HD and with the assumed longitudinal structure of the original color source (which, 
we recall, has support at positive x~ only, with ^ x~ ^ 1/A + ), we shall integrate out the 
quantum fluctuations by using a retarded ie prescription in the LC-gauge propagator (see 
Ref. || for details). We have verified that the same results for \ an d cr would be eventu- 
ally obtained also by using an advanced prescription [^1] . On the other hand, we have not 



been able to give a sense to calculations with principal value or Leibbrandt-Mandelstam 
prescriptions. 

Consider the induced source a first. Up to some tadpoles (i.e., contributions propor- 
tional to / d 2 p±/p 2 \_) which cancel anyway in the final result, this is obtained from 

(W*)> = gf abc ((d + al(x))aUx)) = gTr(T a d+G™(x,y))\ x=y = Tr (T a (Sp(x))) , (13) 

where the equal-time limit of the time-ordered propagator (|T^) is taken as y + = x + + e H. 
Eq. (|T3"|) can be evaluated as follows (more details will be given in Ref. [21]): The relevant 



component of the propagator G n (x, y) (the only one to give a non- vanishing contribution 
to (5p)) can be read off eq. (6.35) of Ref. ||, which gives (omitting tadpoles once again): 

(5p(x)) = -i g ( x \V i -^G ++ V j < i \x), (14) 

where T> % = d l —igA l a T a and T>^ = d^+igA^ with the derivative acting on the function 
on its left. In eq. (|14"D , G ++ (x,y) is a component of the background-field propagator in 
the temporal gauge d~ = 0, which enters at intermediate steps in the construction of the 
propagator in the LC-gauge 0. It is given by eqs. (6.22) and (6.13) of Ref. as: 

G ++ (x,y;p-) = ^Jd 3 zdlG (x-z,p~)5(z-)dlG (z-y,p-) 

x {6(x-)9(-y-)V(x ± )V\z ± ) - 6{-x~)e{y-)V{z x )V\ yi _)} ,(15) 

where G (p) = l/(2p + p~ — p\ + ie). Eq. ( |l~o] ) has been written in the ^"-representation, 
which is convenient given the homogeneity of the problem in time (recall that the back- 
ground fields and sources are independent of x + ; cf. eq. (^)). For the same reason, it is 
preferable to impose the strip restriction (|8|) on p~ rather than on p + [[J : 

A" < |p~| < A~/6, (16) 

where A - = Q\/2A + and Q± is some generic transverse momentum^ 

7 The discrepancy between our results below and those reported in Refs. []l4|, [l^] may be attributed to 
using different gauge-fixing prescriptions. 

8 As explained in M, the two restrictions (H|) and (|l^) are equivalent to leading logarithmic accuracy 
since the relevant quantum contributions involve nearly on-shell quanta, with 2p + p~ m p\ . 
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The expression in eq. (IS) has an intuitive interpretation [| : The particle (a 
"semi-fast" gluon) is moving from y to x in the presence of the singular electric field (0) 
located at z~ = 0. In eq. (|l|), the end points y~ to x~ lie on opposite sides with respect 
to the plane z~ = 0, which corresponds to trajectories crossing the potential^. Consider 
x~ > and y~ < for definitness. On the first part of the trajectory, from y~ to z~ = 0, 
the field A 1 vanishes (cf. eq. flpj)) and the particle moves freely. At z~ = 0, the gluon 
scatters off the electric field (0) with a scattering amplitude V*(zx). Finally, from z" = 
until x~ the color of the gluon is simply rotated, with color matrix V(x±), by the "pure 
gauge" background field A 1 . 

After inserting eq. fllSJ) in eq. (P[) and performing simple manipulations, one obtains: 



M*)) = 9K I d -f- I ifT- e ^ +)x I dr ± ( P± ■ k A _)G ( P )G Q (k) 

J strip 111 p J Z7T Z7T p + + IE J 

{9{jp-)V{x L )V\z L )-e{-p-)V{z L )V\ yA _)}v^ x , (17) 

where it is understood that k~ = p~ (with p~ restricted to the strip (0)), and we have 
used the following shorthand notation for the transverse phase-space integrals: 



dT i = 



(27r) 2 7 (2?r) 



Note the +«e prescription for the "axial" pole at p + = : this is the retarded gauge-fixing 
prescription alluded to before. By using this prescription, together with the Feynman 
prescriptions in the propagators Go(p) and Go(k), the integrals over p + and k + can be 
easily performed via contour technics. This generates ^-functions of ±x~ and ±p~ which, 
together with the functions already present in eq. (|T7j), have the effect to single out 
contributions proportional to 9(x~)9(—p~). 

The final result of this straightforward calculation reads: 

p -ibk + x- _ p-iA+x' f , U 

(5 Pa (x)) = 9 0{x-) Vl / dT ± *^J± k* Tr(T a V z Vj) , (19) 

where the ^"-dependent "form factor" has been generated by integrating over p~ within 
the strip (|16| ) (with k\/2A~ pa A + to the accuracy of interest) : 



i 8, 



dp- 9{-p-) ei J± x - = J_ e - ibA+x ~ - e- iA+x ~ (2Q) 



strip 2n 2(p ) 2 k\ 2irix~ 



Remarkably, eq. ( p!9]) shows that (dp) has support only at positive, and relatively large, 
x~ , with (typically) 1/A + ^ x~ ^ 1/6A + . That is, the induced source is located on 
top on the tree-level source p (with support at x~ ^ 1/A + ). By induction, we conclude 

There are, of course, also "non-crossing" trajectories where x~ and y~ are of the same sign |]. These 
have not been included in eq. (15) since they do not contribute to the induced source (p^). 



7 



that the classical source has support only at x~ > . When the quantum fluctuations 
are integrated out in layers of p + , the source receives new contributions in layers of x~, 
with essentially no overlap between successive layers. This specific longitudinal picture 
is intimately related to our use of the retarded prescription in the LC-gauge propagator. 
With an advanced prescription, one would rather obtain a source with support at negative 
x~ . Still, the integrated (or two-dimensional) charge density o(x±) comes up the same 
with both prescriptions || and reads (cf. eq. fllOD ) 

a s In \ a a (x ± ) = 9 -\ry\v\\ dT ± k* Tr(W^ + ) , (21) 

b ix b J P±k± v ' y±= x J- 

where the covariant derivative involves the asymptotic field: T>\, = d l x — iA 1 ^ (x±). Note 
that the logarithmic enhancement has been generated only after integration over x~ . For 
what follows, it is useful to have the expression of a rotated to the covariant gauge: 

S*i)^i) = 4<% / dvJ-^k^r{T a VlV z ). (22) 

In the weak field limit jaCl, the r.h.s. of this equation scales like gp, as it can be easily 
verified by expanding the Wilson lines; thus a = O(p). In fact, in this limit, a = K, virt p, 
with JC V irt the virtual BFKL kernel [|TJ. In the saturation regime, gat ~ 1, and all the 
non- linear effects displayed by the above equation are equally important. 

Consider \ now. A typical contribution to the charge-charge correlator reads [0, ^ 

(5 Pa (x) 5p b (y)) = 4ig 2 ^(x) G&(s, y) J^ b j (y), (23) 

where it is understood that y + = x + + e. The calculation of this and the remaining 
contributions is similar to that of o outlined above, and eventually gives the following 
result (in matrix notations) 



a 2 if d i d j 
(Sp x 5 Py ) = ^ln-|^(xj^-^|y ± )^ (24) 

+ (p + 2r*(& - (x ± \ ^ \y ± ) (p + 2(vV- ^ V j 



2JJ X ' ^ -Vi V r V 2 

In this equation, both p and T +l are localized at x~ ^ 1/A + . Thus, unlike the induced 
charge ([19|) , the charge-charge correlator appears to be localized near the LC (in both x~ 
and y~)i and thus sensitive to the internal structure of the source. If this was also true 
for Xi h would spoil the separation of scales assumed by the effective theory, and thus 
the validity of the latter. Fortunately, however, this is not the case: Like the electric field 
T +l = d + A l , all the other functions carrying the x~ and y~ dependences in eq. (g4j) are 
total derivatives with respect to x~ or y~ . Thus, x — which is obtained from eq. (|24| ) by 
integrating over x~ and y~, cf. eq. ( |10|) — is sensitive only to the asymptotic fields. 

Specifically, by using the equations of motion ([I]), the following identity can be verified 
(for an arbitrary function <$>(x±)) : 

(p + 2F +i V i ) x $(x ± ) = id+Vm X± ), (25) 
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which allows us to cast x m the following form: 

where all the fields are asymptotic fields, e.g., A l x = ^(xi) and V 1 = d % — igA 1 ^. For 
weak fields, ^4' ~ — <9 J a, and % ~ p)C rea ip, with /C re az the real BFKL kernel After 
rotation to the COV-gauge (cf. eq. (|TT|)), eq. (^) yields: 
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where the derivatives not included in the brackets act on all the functions on their right 
(or left). 

Once x is known, we are in a position to compute also the classical polarization —5a 
in eq. (|Tl"l) . One thus obtains (up to a tadpole which would anyway cancel against the 
tadpole ignored when writing eq. (^)) : 

-gt<Ta(x ± ) = -2idl J dr ± ^-^diTr(T a ViV z ). (28) 

Unlike a in eq. (0), 5a has no contribution linear in p, so it does not contribute at 
the BFKL level. This is as it should, since the BFKL equation is already obtained by 
combining the weak field limits of the previous results for a and x [0 ■ ^ n ^ ne °ther hand, 
5a is as important as a in the non-linear regime where ga ~ 1. By combining eqs. ([22] ) 
and (^), we finally arrive at the following expression for a, which is the quantity which 
enters the RGE for W T [p\ (cf. eq. (0)) : 

ga a (x ± ) = J dT ± ^^Tr(T a VX)j = -gV ± u a (x ± ). (29) 

Up to a normalization factor, the quantity f a (x_i_) introduced above is the induced color 
field in the COV-gauge, that is, the modification in the original field a induced by the 
quantum corrections. It can be verified that u, and therefore a, are free of tadpoles: the 
tadpoles hidden in a, eq. (|2"2"D, and —5a, eq. Q2lf), have cancelled in their sum. 

At this stage, the RGE for W T [p] is fully specified. It has the formal structure shown 
in eq. (H), but with p — > p, x X, an d cr — » a, and x an d & given by eqs. (p7|) and 
(p9|). But the peculiar structure of the "virtual" correction in eq. ( |29| ) suggests that the 
RGE may acquire an even simpler form by using the COV-gauge field a rather than the 
COV-gauge source p as the functional variable of W T . By using (cf. eq. 



f d 2 z ± {x ± \ — \z±) — - — r, (30) 



5p-(x ± ) J -Vi 1 5a a T (z ± ) 

where a^(x±) = a a (x~ = 1/A + , x±) is the color field in the highest bin of rapidity (where 
the quantum corrections are located), one can transform the RGE for W r [p] into the 
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following RGE for W T [a] : 

= a s il l 1 S n d WrVt] - r 7 r [W T <]l , (31) 

dr \2<k^(o; ± )5c^(?/ ± ) L /:cyJ <5a«(x ± ) L xJ J 

where 
and 

g 2 r] ab (x_L,y ± ) = g 2 J d 2 z L J d 2 u L (x±\ \z±) x ab (z±, u±) (u±\ \y ± ) 

The kernels in the equations above have been written in coordinate space by using 
d 2 p± f d 2 k ± p ± ■ k ± Jp± . (x± _ z±) ik± . iz± - y±) 1 (x i - z i )(y i - z l ) 



ip±-(x ± -z ± ) ik ± -(z ± -y ± ) v 'J£ >__ (OA) 

(2tt) 2 J (2tt) 2 p\k\ (2vr) 2 (x ± - z ± ) 2 {y ± - z ± ) 2 ' 1 ; 

The functional derivatives with respect to a T in eq. ( pT|) are understood as variations 
in the color field at x~ ~ 1/A + , which is the end point of the Wilson lines and V (cf. 
the remark after eq. (B)). Thus, the derivatives of the latter read simply: 



8V\x±_ 



I = igT a V\x^) 5 {2 \ Xl _ - z ± ), = -igV( yA _)T^ 2 \y ± - z ± ). (35) 



8a%(z±) ' ~ ~~ ' 5a%(z± 



Eqs. fl3T|)-(|3~3"D represent our main result in this Letter. They govern the flow with 
r = ln(l/x) of the probability density W T [a] for the stochastic color field a T (x) which 
describes the CGC in the COV-gauge. Since the r.h.s. of eq. ( |3"ID is a total derivative 
with respect to a, this flow automatically preserves the correct normalization of the weight 
function: 

J Va W T [a] = 1. (36) 



The functional RGE equation (|3lD can be used to derive evolution equations for all 
the observables which can be related to a (like the soft gluon correlation functions (Q)). 
If 0(a) is any such an operator, then its average over a, defined as in eq. (||) but with p 
replaced by a, obeys the following equation 

Tt (0(o)), = a,(- < ^ {Mv±) + A , (37) 
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which is obtained by multiplying eq. ([Jl|) with 0(a), integrating over a, and performing 
some integrations by parts in the r.h.s. 

In particular, by chosing 0(a) = tr(V^V y ) with the Wilson lines and V y in the 
fundamental representation, and by using eqs. (|32|)-(|35l) and performing the color algebra, 
one eventually obtains the following evolution equation (with N the number of colors) 



d_ 



2vr 2 



d 2 z 



(x± - z±) 2 (y± - z±y 



tr(V}V z ) tr(y2Vy) - Ntr<y}V y )) , (3 



which coincides with the corresponding equation obtained by Balitsky within pertur- 
bative QCD. In the large N limit, the 4-point function in the r.h.s. factorizes: 



tr(VX) tr(Vj V w 



tr(VjV z )) (tr(VX)) (foriV-^oo) 



(39) 



and eq. ( j38l) reduces a closed equation for the quantity Af(x±,y±) = (tr(l — VjV y )) T 
(which represents the forward scattering amplitude of a color dipole off the hadron): 



2^2 



d 2 z 



(x± - y±f 



(x ± - z ± ) 2 (y A 



{n(aT xz + M zy - M xy ) - M xz N zy ] . (40) 



This coincides with the evolution equation obtained by Kovchegov |I7j within the frame- 
work of Mueller's dipole model [18[ . 



We now return to more formal properties of the RGE (|3"T|), and notice the following 
remarkable relation between the coefficients in this equation: 



2 Sr) ab (x ±} y±) 



d z y 



5a b T (y ± ) 



v a (x ± ) 



(41) 



which is most probably a consequence of gauge symmetry (as it relies 1-point and 2-point 
functions). By using this relation, the RGE can be brought into a Hamiltonian form: 



dW T [a] 
dr 



a s 



d 2 x d 2 y 



S 



5a°:(x) 



rj : 



SWr 

Sa b T (y), 



HW T , 



(42) 



with the following, manifestly positive definite, Hamiltonian 

d 2 z ± 



H 

J l a(z±) 



J l a (z ± )J l a (z ± ), 
z % — x % 



2tt 
d 2 x± 

27T (Zj_ - X±) 



id 



x)ab 



5a b T (x ± )' 



(43) 



(Notice that the "current" operator J l a (z±) is indeed Hermitian, since the Wilson lines 
in the adjoint representation are real color matrices: V^ b = Vb a -) This same Hamiltonian 
structure has been recently identified by Weigert [2D| in relation with a generating func- 
tional for Balitsky's evolution equations. Since our effective theory generates Balitsky's 
equations indeed, it should not be a surprise that we come across the same Hamiltonian. 
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Eq. ( f42| ) has the structure of a diffusion equation. Since the Hamiltonian is a purely 
kinetic term, this equation describes a Brownian motion in the functional space of a. 
As r increases, we expect the solution W T [a] to this equation to spread over the whole 
available phase space, and thus to go asymptotically to zero at any given "point" a a (x±) 
(because of the normalization condition (pop). As a simple analogy, consider the equation 
describing the one-dimensional diffusion: 

dW T (x) 1 d 2 W T (x) 

dr 2 V dx 2 ' 1 } 

The solution to this equation, 



.2 



W T (x) = ^e"^, (45) 

is normalized to unity at any r (/ dxW T (x) = 1), and goes smoothly to zero at any x 
when r — > oo, as it spreads over the whole x axis. 



In Ref. p0| ) Weigert has conjectured that a uniform (i.e, a-independent) probability 
distribution W should be a fixed point of the evolution as r — > oo. His arguments are as 
follows: first, any constant Wq is an eigenfunction of H with zero eigenvalue (HWq = 0); 
second, since H is positive definite, all its other eigenvalues must be positive. From 
this, he concludes that the constant eigenfunction (the "fundamental state" of H) will 
be singled out by the evolution for sufficiently large r. However, we disagree with this 
interpretation. As mentioned above, a uniform probability density is necessarily trivial, 
Wo = 0, to cope with the normalization condition. Besides, for the "fundamental state" 
of H to be singled out by the evolution, there must be a gap between this state and the 
first excited state. But such a gap cannot be expected for a Hamiltonian which is just 
a kinetic term, like eq. (fj3|). To distinguish between these various asymptotic scenarios 
and, even more importantly, to find the behaviour of W T [a] at large, but finite r, further 
studies of eq. ([E3) are necessary. 

To conclude, we have established the equivalence between the non-linear evolution 
equations derived within perturbative QCD and within the effective theory for the Color 
Glass Condensate. We have thus resolved a discrepancy between these two approaches 



previously reported in Refs. |14|, [15[ . There are several noticeable differences between our 
calculations and those in Refs. [n, |T3| (chiefly among them, the different gauge-fixing 
prescriptions), which may explain the different results, but a deep understanding of these 
differences requires more study. We have found that the evolution of the gluon correlation 
functions in the non-linear regime at small x is governed by a Hamiltonian describing 
current-current interactions in two spatial dimensions. Given the relative simplicity of 
this Hamiltonian, we believe that analytic solutions to this equation (at least in particular 
limits) should be within our reach. 
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